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ABSTRACT 

Numerical simulations predict that a large number of baryons reside in intergalactic space at tem- 
peratures between 10^ — 10^ K. Highly-ionized metals, such as O VII and O VIII, are good tracers of 
this ^^warm-hot intergalactic medium" , or WHIM. For coUisionally-ionized gas, the ionization fraction of 
each ion peaks at some particular temperature ("peak temperatures"), so different ions can therefore 
trace the IGM at different temperatures. We performed a hydrodynamic simulation to study the metal 
distributions in the IGM. By studying the distribution functions of H- and He-like O, Si and Fe in a 
collisionally-ionized IGM and comparing with semi-analytic results based on the Press-Schechter formal- 
ism, we find: (1) ions with higher peak temperatures (for instance, Fe XXVI) tend to concentrate around 
virialized halos, which can be well described by the Press-Schechter distribution, ions with lower peak 
temperatures are found both in small halos (such as groups of galaxies) and in filaments; (2) lower peak 
temperature ions are more abundant and should be easily observed; (3) peculiar velocities contribute a 
significant part to the broadening of the resonant absorption lines. 



1. INTRODUCTION 

Cosmological N-body simulation plays an important role 
in studying cosmic structure formation and evolution. It 
is of particular interest to study the evolution of bary- 
onic matter because it addresses an important question: 
why is there an apparent decrease in baryon density from 
high to low redshifts (see, e.g., Fukugita, Hogan & Peebles 
1998)? Numerical simulations indicated that the absence 
of baryons at low redshifts represents a yet undetected 
IGM: the "warm-hot intergalactic mediuni^ , or WHIM 
(Dave et al. 2000). For example. Gen & Ostriker (1999) 
and Dave et al. (2000) found that the average temperature 
of baryons is an increasing function of time, with most of 
the baryons at the present time having a temperature in 
the range of 10^ — 10^ K. Although invisible in the optical 
band, the highly ionized metals in WHIM gas might be 
detectable in UV and X-ray bands. 

Since highly-ionized metals, such as O VII and O VIII, 
are good tracers of WHIM gas at 10^ — 10^ K, in this pa- 
per we study the metal distributions in the IGM through 
hydrodynamic simulation. The idea is similar to that 
discussed in Perna & Loeb (1998), Hellsten, Gnedin & 
Miralda-Escude (1998) and Fang & Canizares (2000): 
highly ionized metals would introduce absorption features 
in the X-ray spectrum of a distant quasar — the "X-ray 
Fores f . This paper has the following goals: (1) compare 
predictions of the X-ray forest from semi-analytic models 
to those from numerical simulations, this is important be- 
cause each technique has its own strengths and weakness 
and by contrasting the two, we can determine which re- 
sults are robust; (2) compute predictions for a wider range 
of ions, including O VII, O VIII, Si XIII, Si XIV, Fe XXV, 
Fe XXVI; (3) explore the nature of the gas that gives rise to 
the absorption, including its thermal state and its location 
in clusters, groups or filaments. 
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2. NUMERICAL SIMULATION 

We numerically simulated a cubic region of lOO/i^^ Mpc, 
or 150 Mpc with h = 0.67 on one side at z ^0. We model 
the dark matter with particles (as in a standard N-body 
code) but use the technique of adaptive mesh refinement 
(AMR) to follow the baryons. We refer to Bryan (1996) 
and Norman & Bryan (1999) for further details. In this 
simulation the dark matter particle mass was 5 x 10^^ Mq 
and the smallest cell size (in the high density regions) was 
A9h~^ kpc. To generate spectra, we created a uniform 
grid which contained 512 x 512 x 512 cells, with each cell 
having a size of ~ 195h~^ Kpc, or 19.5 km s~^ in veloc- 
ity space. The simulation itself is a flat low-density cold 
dark matter model (ACDM) with = 0.3, Oa = 0.7 and 
Qb = 0.04. The amplitude of the initial power spectrum 
was fixed so that the rms overdensity in a spherical tophat 
with radius 8h~^ Mpc, as, is 0.9. We use a scale-invariant 
power spectrum with a primordial index of 7i ~ 1. The 
initial cosmological redshift is z = 30. 

Little is known about the metallicity of the IGM. The 
metal abundance could be as low as 10^'^Zq in the low 
column density Lya forest systems (Lu et al. 1998) or 
as high as IZq in the interstellar medium (ISM) of some 
early-type. X-ray luminous galaxies (Matsushita, Ohashi, 
& Makishima 2000). Since we will compare numerical sim- 
ulation with the semi-analytic model in our previous paper 
(Fang & Canizares 2000), we assume uniform metalicity 
and use the metallicities observed from the intracluster 
medium, e.g., 0.5Zq for oxygen, silicon and O.SZq for iron 
(Mushotzky et al. 1996). To investigate X-ray absorption 
in the spectrum of a background source, we need to un- 
derstand the distribution of ionization states of metals in 
the IGM. To obtain the ion density, we assume the gas 
is in coUisional equilibrium and we adopt the ionization 
fractions from Mazzotta et al. (1998). We select helium- 
and hydrogen-like oxygen, silicon and iron. The peak tem- 
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peratures of these ion species range from 10^ K to 10® K, 
which probe different temperature regions of the IGM. 

We have obtained the projected column density, N^, dis- 
tributions of all the six ions. Figure 1 (for O VII) and 
Figure 2 (for Fe XXVI) show two examples. In each figure, 
the ion column density distributions are plotted in four dif- 
ferent panels with > 10^'^, 10^^, lO" and 10^^ cm-^, 
respectively. The color bar on top of each panel indicates 
the logarithm scale from the minimum to the maximum 
column densities in each panel. 
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Fig. 1. — The projected density (column density, iV') distribu- 
tions of O VII, plotted in four different panels with A'^' > 10^'^, 
IQl^, 10" and 10^^ cm , respectively. The color bar on top of 
each panel indicates the logarithm scale from the minimum to the 
maximum column densities in each panel. 

By examining the column density distributions of all 
the six ions, we find that the distributions of these ions 
are closely related to the density and temperature distri- 
butions. Regions with higher baryon densities also have 
higher metal densities. Due to the fact that different 
ions have different peak temperatures, there are differ- 
ences among the distributions of these ion species. For 
instance, many of the filaments that connect high den- 
sity regions are shock-heated to temperatures above 10^ 
K. Since O VII has a peak temperature between 3 x 10"'' 
K and 2 x 10^, we find that a large amount of O VII ex- 
ists in filaments with a column density even higher than 
10^5 (Figure Id). On the other hand, Fe XXVI has a 

peak temperature of over 10^ K. In Figure 2 we find that 
almost no Fe XXVI exists in filaments, and it all concen- 
trates around those collapsed intersections. By comparing 
these two extremes, we find that generally Fe XXVI fol- 
lows the distribution of virialized objects while most of 
O VII exists in filaments, and the distributions of other 
ion species vary between those of Fe XXVI and O VII. 

3. SPECTRAL ANALYSIS 

The methodology we adopt is straightforward: a ran- 
dom line-of-sight (LOS) is drawn across the simulated re- 
gion; temperature, baryon density and velocity are ob- 
tained at each point along the LOS; the ion density is cal- 
culated based on the metallicity and ionization fraction; 
the synthesized spectrum is then obtained from the opti- 



cal depth along the LOS; finally, by fitting the resonant 
absorption lines in the simulated' spectrum we obtain 
ion column densities (iV*) and the Doppler parameters {b). 
This method is similar to the method used in the Lya for- 
est; we refer to Zhang et al. (1997) for a complete discus- 
sion and formalism. We draw a total of 5000 random LOS 
across the simulated region. For each LOS we generate 
a simulated spectrum and fit it by the "threshold-based' 
method. We then calculate the so-called "X-ray forest" 
distribution fmictioii (XFDF), d'^P/dN^dz, defined as the 
number of absorption systems along the line of sight with 
a column density between and + dN^ per unit red- 
shift. 





N.{Fe XXV) > 10" cm-' N,(rB XXV) > 10" em"° 





H,(Fe XXV) > 10" cm-" N,{re XXV) > 10" cm"* 

Fig. 2. — The projected column density distributions of Fe XXVI. 

Based on the formalism in Fang & Canizares (2000), 
we compare our numerical simulation with the analytic 
model. The Press-Schechter function (Press & Schechter 
1974) predicts the spatial distribution of virialized objects. 
We would expect that the analytic model and the numeri- 
cal simulation should give consistent results for the distri- 
bution of the high-peak-temperaturc ions, such as Fe XXV 
and Fe XXVI, because these ions follow the distribution of 
virialized objects; on the other hand, we would not expect 
the low-peak-temperature ions, such as O VII and O VIII, 
would have the same behavior because most of these ions 
are distributed in filaments. 

Figure 3 shows the XFDF predicted from the Press- 
Schechter approach for all the six ions (solid line in each 
plot) and also displays the distributions from the numer- 
ical simulation (the fiUed-dottcd line). The most striking 
feature is that the predicted distributions fit the simu- 
lation very well at high column density regions, specif- 
ically at Ni > 10^^ cm~^ for both oxygen species and 
at Ni > 10^''' cm^^ for silicon species and Fe XXV. For 
Fe XXVI, both distributions agree very well in the whole 
column density range (10^^ — 10^^ cm~^). This result can 
be well understood in the framework we describe at the be- 
ginning of this section: high-peak-temperature ions, such 
as Fe XXVI, follow the distribution of collapsed, virialized 
halos well described by the PS formalism. For low-peak- 
temperature ions, the high column densities also come 
from virialized structures, while those with low column 
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density are distributed in those filaments which can not 
be described by the PS distribution. Thus except for 
Fe XXVI, at low column density the PS distributions are 
systematically lower than in the numerical simulation. 




(cm ') 



N, (cm-') 



10l2 10l3 10l4 10l5 10I6 lol 
(cm ") 




N, (cm-") 




Ni (cm"" 



N, (cm-") 



Fig. 3. — The XFDF predicted from the Press-Schechter approach 
(solid lines) and numerical simulation (fiUed-dotted lines). 

Another feature seen in this figure is that all curves 
display a sharp cutoff somewhere between fO^^ and 
10^'' cm~^. The reason for this apparent drop at high 
column density end is two fold: first, systems with high 
column density generally correspond to high temperature 
regions, and at high temperature the ionization fraction 
drops quickly and so does the ion column density; second, 
high colmnn density systems are closely related to virial- 
ized objects, such as groups and clusters of galaxies, and 
the Press-Schechter theory predicts that the spatial den- 
sity of virializcd objects decreases exponentially as their 
masses (and temperatures) increase. 

Doppler 6-parameter governs the line width and can be 
affected by several factors, such as gas temperature (ther- 
mal motion), the Hubble flow and the velocity structure. 
Figure 4 shows the 6-paramctcr distribution function for 
six ions (the filled circles). All six ions exhibit a low cutoff 
to the Doppler parameter. The cutoffs are primarily due 
to the cutoff temperatures for each ionization. 

To understand the role of different line broadening 
mechanisms, we first fit the ^-parameter distribution with 
a Gaussian. We find that there exists a non-Gaussian tail 
in each plot. This tail can be well fitted by a lognormal 
distribution (the solid line): / oc exp [—log (fe/6;„)/2(T;^„] 
Here bin and are the mean velocity and the standard 
deviation of this lognormal distribution. We find this dis- 
tribution naturally cuts off at low b, and also provide an 



excellent fit to the long tail. This non-Gaussian tail im- 
plies that the Hubble flow and peculiar velocity contribute 
a significant part to those ions with high b values. 
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Fig. 4. — The &-paj:ameter distribution function for six ions (the 
dotted lines), fitted with Gaussian distributions (the dashed lines) 
and log normal distributions (the solid lines). 

4. SUMMARY AND DISCUSSION 

In this paper we have performed a hydrodynamic sim- 
ulation based on a ACDM model. By drawing random 
LOS and synthesizing the X-ray spectrum for each LOS, 
we obtained the XFDFs for hydrogen and helium-like O, 
Si, Fc. We have investigated the column density distribu- 
tions of highly-ionized metals by comparing the numeri- 
cal simulation to the analytic method described in Fang 
& Canizares (2000). We also discuss the distribution of 
Doppler 6-parameters. Our main conclusions can be sum- 
marized as follows: 

1. Ions with low peak temperatures such as O VII 
and O VIII trace small halos (galaxy groups) and 
filaments, while high peak ions such as Fe XXV and 
Fe XXVI are only observed in the largest objects 
(galaxy clusters). For all ions, the strongest signal 
come from virialized halos; however, filaments 
show up as lower column density absorbers. 

2. The low peak temperature ions are more abundant 
than ions with high peak temperatures. This can 
be understood from numerical simulations that 
most of the baryons have temperatures ranging 
from 10^ — 10^ K and the baryon fraction decreases 
as the temperature increases. This implies that 
from the point view of observation, ions with 
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low peak temperatures such as O VII and O VIII 
should be more easily observed. 

3. We compared the column density distribution of 

absorbers predicted by the simulations against an 
analytic model developed in Fang & Canizares 
(2000) and found that the numerical simulation 
confirms our predictions: ions with high peak 
temperatures (T > lO^K) follow the distributions 
of collapsed halos, which can be described by the 
Press-Schechter distribution; for ions with low 
peak temperatures, those lines with high column 
density come from virialized structures, while 
those with low column density are distributed in 
diffuse filaments which are not well described by 
the PS distribution. 

4. Higher excitation ions produce lines which are 
generally wider (as determined by the Doppler 
6-parameter) than low excitation lines, although 
there is a wide distribution in all cases. This 
distribution is well fit by a log normal distribution. 
The sharp lower cutoff of this distribution is 
provided by the fact that there is a minimum 
temperature required to create the ion in the first 
place. Peculiar velocities in the gas also contribute 
significantly, providing a non-Gaussian tail at the 



high-6 end. 

The largest uncertainty in spectral synthesis comes from 
the metallicity: we did not include any gas processes in our 
numerical simulation. Although these gas processes are 
not important in heating and cooling gas with T > 10^ K, 
the IGM metals enriched by these processes are crucial in 
UV/X-ray observations: it is these metals which produces 
detectable features. We assume a uniform metallicity in 
generating the synthesized spectrum, however, we would 
expect that metallicity should be determined by the metal 
enrichment history of the IGM, which in turn depends on 
a variety of sources: the UV/X-ray background radiation, 
supernova input, etc, and these are likely to be different in 
low and high density regions. It is important to incorpo- 
rate all these processes to obtain a self-consistent result of 
metal distributions and ionization structures in the IGM. 
This is the goal of our future numerical simulations. 
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